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ABSTRACT 



Astronomical Institute of the Academy of Sciences of the Czech RepubHc, v.v.i. 
Fricova 298, Ondrejov, CZ-251 65 



2. INFORMATION OBTAINED 
TROSCOPIC ANALYSIS 



FROM SPEC- 



Despite the lack of important functions supported only 
by legacy applications, the current VO-compatible tools 
have enough capabilities to allow powerful analysis of 
stellar spectra using both public archives and local pro- 
prietary data. We give examples of the possible multi- 
wavelength analysis of changes of profiles in different 
spectral lines using the current SSAP-compatible VO 
tools. The potential of future VO applications supporting 
more elaborate methods for stellar analysis is discussed 
as well. 



Key words: spectroscopy, spectral lines, Doppler Imag- 
ing, SSAP , Virtual Observatory. 



1. INTRODUCTION 



Astronomical spectroscopy uses a wide range of tech- 
niques with different level of complexity to achieve its 
final goal — to estimate the most precise and rehable in- 
formation about celestial objects. 

The scientific analysis of such spectra requires further 
processing by the variety of different methods. In cer- 
tain studies a huge number of spectra has to be collected 
from different servers (in different spectral regions) and 
transformed into common units. 

Accomplishing the multi-spectral analysis in VO en- 
vironment may benefit from automatic aggregation of 
distributed archive resources, seamless on-the-fly data 
conversion, common interoperability of all tools (using 
PLASTIC or SAMP protocol) and powerful graphical vi- 
sualisation of measured and derived quantities 



The spectral lines give us a plenty of information about 
the physical nature and evolution of the given astronom- 
ical object. Examples of information estimated from dif- 
ferent parameters of spectral lines are given below: 

• Position of line (wavelength) 

- Individual chemical elements present in the 
stellar atmosphere 

- Excitation / ionisation state of given element 

- Structure of molecules (e.g. rotational and vi- 
brational states of molecules in IR regions) 

- Radial velocity (RV). If variable RV is mea- 
sured, it may be the sign of binarity, some- 
times even the estimation of orbital parameters 
is feasible 

• Line Shape 

- Stellar parameters (Toft, log g, rotation) 

- Stellar activity (turbulence, granulation) 

- The profiles of core/wings can present differ- 
ent physics dependent on optical depth (e.g. 
limb darkening law) 

- Expansion of gas shells, winds (P Cyg profiles 
in Novae, double-peaked emission profiles in 
Be stars) 

• Line profile variability in time 

- Change of physical state (e.g trigger of of 
emission phase in Be stars, outbursts of novae) 

- Stellar spots (magnetic field structures, local 
overabundance - typical for Ap stars) 

- Pulsations (S Cep, RR Lyr, Miras) 

- Non radial pulsations (NRP) (e.g. 5 Set, /? Cep 
types) 

- Detection of extra-solar planets in spectra (bi- 
sector method makes small contributions in 
global line profile enlarged and variable during 
planet's transit) 
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3. SPECTRA POST-PROCESSING 



The current VO applications using SSA protocol are very 
simple providing only the URL of the given dataset, so 
the client gets a whole (often quite large) file. The fur- 
ther processing is then done fully by the client which 
has to download all the relevant datasets and keep them 
in the memory or local storage. The real strength of 
the VO technology lies, however, in transferring part 
of the client's work to the server, which usually runs 
on a powerful machine with fast connection to the data 
archives. The typical examples of commonly required 
post-processing of fully reduced spectra (at least in stel- 
lar astronomy) are given below: 

• Cutout services (selection of only certain spectral 
lines or regions within the given wavelength range) 

• Projection of multidimensional datasets (in 3D spec- 
troscopy) 

• Rectification (normalisation) of continuum 

• Rebinning to given, usually equidistant grid of 
wavelengths (constant AA or A In A) 

• (De)convolution of instrumental profile 

• Application of physical broadening functions (rota- 
tion, limb darkening) 

• Shift in radial velocity, application of heliocentric 
correction computed on server 

• Merging of individual echelle orders in a single 
spectrum 



Although most of these post-processing methods can be 
implemented in a straightforward way, the rectification of 
continua may be a difficult problem, especially on echelle 
spectra. 
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Figure 1. Dynamical spectrum of H/3 line profile vari- 
ability of 59 Cyg. Residuals from average profile of 38 
spectra are binned to 20 phase bins corresponding to pe- 
riod 28.192 days. Two expanded cycles shown for clarity. 
Individ ual profiles are overplotted above. After \MaintA 



averaged one. First the average of many high disper- 
sion high SNR spectra (with removal of outliers) is pre- 
pared (called template spectrum). Then each individual 
spectrum in time series is either divided by the template 
(quotient spectrum) or the template is subtracted from 
it (the differential spectrum). The group of similar re- 
sulting intensities is given the same colour or level of 

gr ay. See Fig. 11 More examples may be found e.g. 

inlde Jong et 31.1119991) . iMaintz. (,2003h or lUvtterhoevenI 
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4. ADVANCED SPECTRAL ANALYSIS 



As the goal of the VO is to make easy and comfortable the 
physical analysis of a huge number of fully reduced (and 
post-processed) spectra in the environment of VO client 
(or web portal), the common recipes of spectral analy- 
sis have to be implemented as VO-compatible. The ba- 
sic techniques were reviewed in ISkodal (l2008l) . we will 
briefly mention only the more complicated tasks of spec- 
tra analysis: 



4.2. Measurement of radial velocity and higher mo- 
ments of Une profile 

The one of the important information received from spec- 
trum is the radial velocity (RV) of the object. From its 
changes the binarity can be revealed, or the possession of 
extrasolar planet. The combination of higher moments of 
line profile is a one of the possible ways of de termination 
ofnon radial pulsation modes - numbers I , m dAerts et al.l 
119921) . 



4.1. Dynamical spectrum 



4.3. Measurement of equivalent width 



It is sometimes called the gray representation or trailed 
spectrum. The basic idea is to find the small time- 
dependent deviations of individual line profiles from an 



The Equivalent width (EW) of the spectral line gives the 
information about the number of absorbing or emitting 
atoms of given element. Emission lines by definition have 
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negative EW. The changes in EW during time may bring 
about the information about the dynamic evolution of the 
target and may be subjected to period analysis. 



4.4. Bisector analysis 



It is a method describing quantitatively the tiny asymme- 
try or subtle changes in line profiles. The characteris- 
tic shape of bisector gives the information about turbu- 
lence fields (e.g. convection) in stellar photosphere, char- 
acterised b y the value of micro-turbulent velocity (iGravl 
Il982ll2005i) or about other processes causing the tiny pro- 
file asymmetry. It has be en used successfull y for search- 
ing of extrasolar planets (iPovich et al.ll200ll) or in astero- 
seismology. It requires high resolution normalised spec- 
tra with extremely high SNR. 



4.5. Period Analysis 



Its aim is to find the hidden periods of variability of given 
object. Sometimes this period can be identified with some 
physical mechanism (e.g. orbital period of binaries, rota- 
tional modulation or pulsations). Wide range of objects 
show the multi-periodicity on various time scales (e.g. bi- 
nary with pulsating components). Ones the suspected pe- 
riod is found, the data may be folded accordingly, plotted 
in circular phase corresponding to this period. 



4.6. Doppler imaging 



It was introduced by 'Vogt & Penrod' ('1983a') as a method 
allowing the surface mapping of stellar spots. First test 
were done on stars o f RS CVn type and on Oph 
dVogt & Penrodll983bl) . Works well on rapid rotators and 
needs a high resolution spectra with very high SNR (300- 
500). Tho whole rotational period should be covered 
well, better several times. When all the requirements are 
met, the map of surface features (spots, nodes of non ra- 
dial pulsations) is obtained with very high accuracy. See 
the left panel of Fig.|2l 




Figure 2. Zeema n Doppler Imaging of II Peg. After 
\StrassmeieA f l2Q0g|) 



4.8. Spectra disentangling 

This method allows to separate the spectra of individual 
stars in binary or multiple systems and simultaneously 
to find orbital parameters of the system, even in case of 
heavy blending of lines. It supposes the changes in line 
profile are caused only by combination of Doppler shifted 
components (no intrinsic variability of star). The best so- 
lution of orbital parameters and disentangled line profiles 
of individual stellar components are found by least square 
global minimisation. The method also enables the re- 
moval of the telluric lines with great precision. The good 
orbital coverage and the estimate of orbital parameters is 
required. Two approaches exist: 

Wavelength s pace disentangling 

developed bylBagnuolo & GiesI (Il991h and improved by 
ISimon & Sturm' ('1994'). It needs a large memory to 
store sparse matrices, requires large computing power. 
It is more straightforward to understand the results and 
sources of errors. 



Fourier space disentangling 

introduced by Hadrava ( 1995l ll997h in program KOREL. 
Another program avail able today (stil l based on KOREL 
ideas) is FDBINARY dlliiic et al .112004) . They work in 
Fourier space, and transform the wavelengths into In A. 
They solve a small amount of linear equations, so they are 
memory savvy and can be run on even small computer. 
The method, however, requires perfect continuum fit. 



4.7. Zeeman Doppler Imaging 



4.9. Classification of stellar spectra 



Quite complicated processing of spectra is required for 
study of stellar magnetic fields. The estimation of 
magnetic field from polarimetry using the Zeeman phe- 
nomenon involves the processing of long series of homo- 
geneous spectra to be accomplished in parallel with ex- 
treme precision and requires again the information from 
synthetic models (simulation of Stokes parameters on 
simulated m agnetic stars) The n ice example is the model 
of II Peg bv IStrassmeieij (l2008h . See the right panel of 
Fig.El 



It used to be a very time-consuming and quite subjec- 
tive method in the era of photographic plates (e.g. MKK 
classification). Today it has been done finding the mini- 
mal differences between a grid of template (or even syn- 
thetic) spectra and the examined one. The techniques of 
artificial intelligence (e.g. neural networks) may be very 
efficient, but the global optimisation using the min- 
imisation seems to be more reliable with possibility of 
iterative control. An example of automati c classificat ion 
engine for white dwarfs is described bv i Winter! (l2006h . 
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5. KILLER VO SPECTRAL APPLICATIONS 



Here we give some scientific cases that may become the 
killer application forcing the wide astronomical commu- 
nity to use the VO tools with capabilities suggested above 
despite its current reluctance and hesitation: 

• Use VO to find all stars with emission in given line, 
find the time when it was in emission and plot the 
time evolution of its EW. 

• Use VO to get 1000 or more spectra of the given 
object, cut out regions around given lines, plot the 
lines, make a gray dynamic spectrum folded in time 

• The same, but fold by period clicked on interactive 
periodogram 

• Get the unknown lines identification of piece of 
spectra from theoretical observatory (SLAP - Sim- 
ple line access protocol) having the line selection 
limited by pre-estimated temperature (using Saha 
equation) 

• Create light and radial velocity curve of a binary star 
for given period (estimated by other VO tool running 
in parallel and exchanging data over PLASTIC or 
SAMP) 

• Fit the grid of models (Toff, log g) to the observed 
spectrum for many stars. 

• Extract detailed profiles of given spectral lines from 
archives of echelle spectra where every echelle or- 
der is kept separately (unmerged spectra in so called 
pixel-order space) 



6. CONCLUSIONS 



The large part of spectroscopic analysis today has 
been accomplished by several independent non VO- 
compatible legacy packages as the currently available 
VO clients supporting the SSA protocol can provide with 
only simple, mostly interactive capabilities on individual 
spectra. 

The processing of large amount of spectra in a consistent 
way requires the simple post-processing (e.g. cutout or 
rebinning) services to be implemented on the server side 

A lot of interesting physical information about astronom- 
ical objected can be obtained from only several short 
spectral ranges covering selected spectral lines. Dynamic 
behaviour of interesting objects can be estimated from the 
changes of line profiles as well as from period analysis of 
their moments (including equivalent width). 

By introduction of modern VO-aware tools into the astro- 
nomical spectral analysis a remarkable increase of effec- 
tiveness of astronomical research can be achieved. 
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